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Abstract 

Scalarization  approaches  to  purposely  generating  epsilon-efficient  solutions  of  multiobjective 
programs  are  investigated  and  a  generic  procedure  for  computing  these  solutions  is  proposed 
and  illustrated  with  an  example.  Real-life  decision  making  situations  in  which  the  solutions  are 
of  significance  are  described. 

Keywords:  multiobjective  programs,  e-efficient  solutions,  e-Pareto  outcomes,  e-nondominated 
outcomes,  scalarizations,  e-optimality. 


1  Introduction 

In  the  middle  of  the  nineteen  eighties,  Loridan  (1984)  introduced  a  notion  of  e-efficient  solutions 
for  multiobjective  programs  (MOPs),  which  was  followed  by  White  (1986)  who  proposed  several 
concepts  of  approximate  solutions  for  MOPs  and  drafted  methods  for  their  generation. 

For  the  last  two  decades,  e-efficient  or  approximate  solutions  of  MOPs  have  been  examined 
in  the  literature  by  many  authors  from  different  points  of  view.  Existence  conditions  were  devel¬ 
oped  by  Deng  (1997)  and  Dutta  and  Vetrivel  (2001)  for  convex  MOPs  while  KKT-type  conditions 
and  saddle  point  conditions  were  derived  by  Dutta  and  Vetrivel  (2001)  and  Liu  (1996).  The  lat¬ 
ter  also  proposed  e-properly-efficient  solutions  for  convex  non-differentiable  MOPs.  Connections 
between  different  definitions  of  approximate  solutions  were  analyzed  by  Yokoyama  (1996,  1999). 
Valyi  (1985),  Tarnmer  (1994),  Tanaka  (1996),  and  others  studied  approximate  solutions  of  vector 
optimization  problems  in  general  ordered  vector  spaces. 

Approximate  efficient  solutions  have  also  been  used  for  approximation  of  the  Pareto  set  of  spe¬ 
cially  structured  biobjective  programs  (BOPs).  In  this  vein,  Ruhe  and  Fruhwirth  (1990)  developed 
an  approximation  algorithm  for  a  minimum  cost  flow  problem,  Blanquero  and  Carrizosa  (2002) 

‘Partially  supported  by  the  Automotive  Research  Center,  a  U.S.  Army  TACOM  Center  of  Excellence  for  Modeling 
and  Simulation  of  Ground  Vehicles  at  the  University  of  Michigan. 
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worked  with  a  location  problem,  and  Angel  et  al.  (2003)  studied  a  scheduling  problem.  Further¬ 
more,  White  (1998b)  applied  the  concept  of  e-efficient  solutions  to  a  portfolio  selection  problem, 
and  to  a  linear  MOP  with  a  specially  structured  feasible  set  (White  (1998a)).  Fadel  et  al.  (2002) 
attempted  to  examine  the  curvature  of  the  Pareto  curve  for  BOPs  using  e-efficiency  as  a  measure 
of  sensitivity. 

In  view  of  the  literature,  the  current  belief  is  that  the  concept  of  e-efficient  solutions  accounts 
for  modeling  limitations  or  computational  inaccuracies,  and  thus  is  tolerable  rather  than  desirable. 
Consequently,  methods  purposely  avoiding  efficiency  and  guaranteeing  e-efficiency  have  not  been 
well  developed. 

The  objective  of  this  paper  is  to  propose  methods  for  generation  of  e-efficient  solutions.  In 
Section  2,  we  formulate  the  problem  of  interest  and  present  the  terminology  used  in  this  paper. 
In  Section  3,  we  explore  relationships  between  MOPs  and  associated  single  objective  programs 
(SOPs),  a  subject  matter  that  has  earlier  been  undertaken  by  some  authors.  References  specific 
for  each  of  the  SOPs  are  listed  in  the  subsections  of  Section  3  in  which  these  SOPs  are  examined. 

White  (1986,  1998b)  seems  to  be  the  only  one  to  have  addressed  the  practical  issue  of  computing 
e-efficient  solutions.  The  scheme  he  advocates  consists  in  the  computation  of  approximate  efficient 
solutions  by  means  of  approximate  optimal  solutions  of  an  SOP  associated  with  the  original  MOP. 
In  Section  4,  we  follow  on  this  scheme  and  offer  specific  guidance  relating  approximate  optimality 
to  approximate  efficiency.  In  Section  5,  we  present  an  illustrative  example  and  discuss  real-life 
decision  making  situations  that  motivate  this  study  of  approximate  efficiency.  Section  6  concludes 
the  paper. 

2  Terminology  and  Problem  Formulation 

We  first  introduce  some  basic  notation.  Let  Mn  and  Rm  be  finite  dimensional  Euclidean  vector 
spaces  and  y,y'  €  Mm.  y  >  y'  denotes  y,;  >  y'i  for  all  i  =  1, . . . ,  m.  y  ^  y'  denotes  r/j  >  y[  for 
all  i  =  1, . . . ,  m.  y  >  y'  denotes  y  ^  y'  but  y  /  y' .  The  relations  <  and  <  are  defined  in  the 
obvious  way.  Let  :=  {y  G  Rm  :  y  ^  0}.  The  sets  are  defined  accordingly. 

Let  X  C  Rn  be  a  feasible  set  of  decisions  and  let  /  be  a  vector-valued  objective  function 
/  :  Mn  — *  Wn  composed  of  m  real- valued  functions,  /  =  (/i,...,/m),  where  /)  :  M"  — ■>  M  for 
i  =  1, . . . ,  m.  The  multiobjective  program  (MOP)  is  given  by 

MOP:  minimize  (/i(x), . . . ,  fm(x))  subject  tor£l 
where  the  minimization  is  understood  as  finding  the  set  of  efficient  solutions  in  X. 

Definition  2.1  Consider  the  MOP.  A  point  x  £  X  is  called 

(i)  a  weakly  efficient  solution  if  there  does  not  exist  such  that  f(x)  <  f(x); 
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(ii)  an  efficient  solution  if  there  does  not  exist  x  E  X  such  that  f(x)  <  f(x). 

A  feasible  solution  (decision)  x  E  X  is  evaluated  by  the  m  objective  functions  producing  the 
outcome  f(x).  We  define  the  set  of  all  attainable  outcomes  or  criterion  vectors  for  all  feasible 
solutions  in  the  objective  space,  Y  :=  f(X)  C  Rm.  The  image  f(x)  €  Y  of  a  (weakly)  efficient 
point  is  called  a  (weak)  Pareto  outcome. 

Yu  (1974)  introduced  the  concept  of  nondominated  outcomes  as  a  generalization  of  Pareto 
outcomes. 

Definition  2.2  Let  D  C  Mm  be  a  cone  and  Y  C  Mm.  Then  y  EY  is  called 

(i)  a  nondominated  outcome  of  the  MOP  if  there  does  not  exist  y  E  Y,y  y  and  d  E  D  such  that 
y  =  y  +  d,  or  equivalently,  (y  —  D\  {0})  fl  Y  =  0; 

(ii)  a  weakly  nondominated  outcome  of  the  MOP  if  (y  —  int  D)  n  Y  =  0. 

For  D  =  M>,  the  nondominated  outcomes  become  the  Pareto  outcomes.  Following  Loridan 
(1984)  we  define  e-efficient  solutions  of  the  MOP. 

Definition  2.3  Consider  the  MOP  and  let  e  G  .  A  point  x  G  X  is  called 

(i)  a  weakly  e-efficient  solution  if  there  does  not  exist  x  G  X  such  that  f(x)  <  f(x)  —  e; 

(ii)  an  e-efficient  solution  if  there  does  not  exist  x  G  X  such  that  f(x)  <  f(x)  —  e. 

The  image  f(x)  E  Y  of  a  (weakly)  e-efficient  point  is  called  a  (weak)  e-Pareto  outcome.  Fol¬ 
lowing  Yu  (1974)  we  define  (weakly)  e-nondominated  outcomes. 

Definition  2.4  Let  D  C  Rm  be  a  cone,  e  E  D,  and  Y  C  Mm.  Then  y  EY  is  called 

(i)  an  e-nondominated  outcome  of  the  MOP  if  there  does  not  exist  y  E  Y,y  /  y  and  d  E  D  such 
that  y  =  y  +  d  +  e,  or  equivalently,  (y  —  (D  +  e)  \  {0})  0  Y  =  0; 

(ii)  a  weakly  e-nondominated  outcome  of  the  MOP  if  (y  —  int(D  +  e))  O  Y  =  0. 

Given  the  MOP,  we  formulate  a  scalarized  (single  objective)  program  (SOP).  Let  S  C  X  be 
a  subset  of  the  decision  set  X,  U  be  a  set  of  auxiliary  variables,  and  II  be  a  set  of  parameters. 
Let  T  :=  f(S)  denote  the  set  of  attainable  outcomes  for  the  SOP  and  s  :  T  x  U  x  II  — >  3R  be  a 
scalarizing  function.  Then  the  SOP  associated  with  the  MOP  is  given  by 

SOP(7r):  minimize  s(f(x),u,  ir)  subject  to  x  E  S,u  E  U 

where  ir  E  II  is  a  vector  of  parameters  chosen  by  a  decision  maker.  We  next  review  the  notion  of 
optimality  and  e-optimality  for  the  SOP. 
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Definition  2.5  Consider  the  SOP  and  let  e  >  0.  A  point  (x,u)  G  S  x  [/  is  called 

(i)  an  optimal  solution  if  s(f(x),u,n)  <  s(f(x),u,  ir)  for  all  (x,u)  G  S  xU;  the  outcome  y  = 
f(x)  G  T  is  called  optimal; 

(ii)  strictly  e- optimal  solution  if  s(f(x),u,n)  <  s(f(x),u,  7r)+e  for  all  (x,u)  G  SxU;  the  outcome 
y  =  f(x)  G  T  is  called  strictly  e-optimal; 

(in)  e-optimal  solution  if  s(f(x),u,n)  <  s(f(x),u,n)  +  e  for  all  (x,u)  G  S  x  U;  the  outcome 
y  =  f(x)  G  T  is  called  e-optimal. 

Note  that  by  definition,  each  optimal  solution  of  the  SOP  is  in  particular  strictly  e-optimal,  and 
each  strictly  e-optimal  solution  is  also  e-optimal. 

3  e-Optimality  and  e-Efficiency 

In  this  section,  we  review  seven  well-known  scalarization  methods  for  the  MOP  and  establish  the 
corresponding  relationships  between  e  G  M  and  e  G  Mm  when  solving  both  the  scalarized  problem 
and  the  MOP  for  e-optinrality  and  e-efficiency,  respectively. 

3.1  Weighted-sum  scalarization 

The  weighted-sum  method  for  multiobjective  programming  is  defined  as 

m 

WS(ro):  minimize  Wjfj(x)  subject  to  x  G  X, 

1=1 

where  w  G  R>  is  a  given  weighting  parameter.  Observe  that  for  this  scalarization,  the  feasible  set 
is  given  by  S  =  X  so  that  T  =  Y  is  the  complete  set  of  outcomes.  Moreover,  U  =  0  and  II  =  . 

White  (1986),  Deng  (1997)  and  Dutta  and  Vetrivel  (2001)  studied  this  method  in  the  context 
of  relationships  between  e-optinrality  for  the  SOP  and  e-efficiency  for  the  MOP.  Following  on  their 
work,  we  develop  more  specific  results. 

Proposition  3.1  Given  the  MOP  and  e  G  Rm,  let  e  <  E?;=i  wi£i- 

(i)  Ify&Y  is  strictly  e-optimal  for  the  WS(w),  then  y  is  e-Pareto  for  the  MOP. 

(ii)  Ify&Y  is  e-optimal  for  the  WS(w),  then  y  is  weak  e-Pareto  for  the  MOP. 

(in)  If  y  €  Y  is  e-optimal  for  the  WS(w)  with  w  G  then  y  is  e-Pareto  for  the  MOP. 

Proof  Although  quite  similar,  we  give  all  three  proofs  in  detail. 
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For  (i),  let  y  G  Y  be  strictly  e-optimal  for  the  WS(u>), 

m  m 

y  mm  <  y  WiVi  +  e  for  a11  y  e  v  / 

2—1  2=1 

and  suppose  that  y  is  not  e-Pareto.  Then  there  exists  y  G  Y,y  ^  y  such  that  y  >  y  +  e,  yielding 

m  m  m  m 

y  Wiyi  >  y  mm + y  wi£i  >  y  mm + e 
2=1  2=1  2=1  2=1 

in  contradiction  to  the  above.  Hence,  y  is  e-Pareto  for  the  MOP. 

For  (ii),  let  y  G  Y  be  e-optimal  for  the  WS(iu), 

m  m 

y  mm  <  y  myi  +  e  for  all  y  G  Y,  y  /  y, 

2=1  2=1 

and  suppose  that  y  is  not  weak  e-Pareto.  Then  there  exists  y  £  Y,y  /  y  such  that  y  >  y  +  e, 
yielding 

m  m  m  m 

y  mm  >  y  my%  +  y  Wi£i  >  y  wm  +  e 
2=1  2=1  2=1  2=1 

in  contradiction  to  the  above.  Hence,  y  is  weak  e-Pareto  for  the  MOP. 


For  (iii),  let  y  G  Y  be  e-optimal  for  the  WS(w)  with  w  G  M>, 

m  m 

y  mm  <  y  mVi  +  e  for  all  y  G  Y,  y  ^  y, 

2=1  2=1 

and  suppose  that  y  is  not  e-Pareto.  Then  there  exists  y  G  Y,  y  /  y  such  that  y  >  y  +  e,  yielding 

m  m  m  m 

y  wm.  >  y  Wiyi  +  y  Wi£i  >  y  Wiyi  +  e 

2=1  2=1  2=1  2=1 

in  contradiction  to  the  above.  Hence,  y  is  e-Pareto  for  the  MOP.  □ 

Corollary  3.1  Given  the  MOP  and  e  G  Rm,  let  e  <  wi£i- 
(%)  If  x  G  X  is  strictly  e-optimal  for  the  WS(w),  then  x  is  e-efficient  for  the  MOP. 

(ii)  If  x  G  X  is  e-optimal  for  the  WS(w),  then  x  is  weakly  e- efficient  for  the  MOP. 

(iii)  If  x  G  X  is  e-optimal  for  the  WS(w)  with  w  G  R>,  then  x  is  e-efficient  for  the  MOP. 
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3.2  Constrained-objective  scalarization 

The  constrained-objective  scalarization  of  the  MOP,  also  known  as  the  epsilon-constraint  method 
(see  Chankong  and  Hainres  (1983)),  is  defined  as 

CO((5):  minimize  fi(x)  subject  to  x  £  S  :=  {x  £  X  :  fi(x)  <  <5j, i  =  2, . . . ,  m}, 

where  <5*  £  M,  i  =  2, . . . ,  m  are  given  upper  bounds  on  the  objectives  ft,  i  =  2, . . . ,  m.  Observe  that 
for  this  scalarization  T  =  {y  £  Y  :  yi  <  5i,  i  =  2, . . . ,  m U  =  0  and  II  =  Rm_1. 

Relationships  between  e-optinral  solutions  of  this  scalarization  and  e-efficient  decisions  of  the 
MOP  were  also  examined  by  White  (1986). 

Proposition  3.2  Given  the  MOP  and  £  £  let  e  <  £\. 

(i)  If  y  £  T  is  strictly  e- optimal  for  the  CO(5),  then  y  is  e-Pareto  for  the  MOP. 

(ii)  If  y  £  T  is  e-optimal  for  the  CO(S),  then  y  is  weak  e-Pareto  for  the  MOP. 

Proof  We  only  prove  (i)  and  then  (ii)  follows  analogously.  Let  y  £  T  be  strictly  e-optimal  for  the 

CO  (d), 

111  <  yi  +  e  for  all  y  &T.y^  y, 

and  suppose  that  y  is  not  e-Pareto.  Then  there  exists  y  £  Y,  y  /  y  such  that  y  >  y  +  e,  yielding 

yi>yi  +  £i>yi  +  e; 

in  >  yi  +  £i  >  yi,  i  =  2, . . . ,  m. 

Since  y  £  T,  we  have  fji  <  Si,  i  =  2, . . . ,  m,  and  therefore  y  £  T  in  contradiction  to  the  above. 
Hence,  y  is  e-Pareto  for  the  MOP.  □ 

Corollary  3.2  Given  the  MOP  and  e  £  M>  ,  let  e  <  £\. 

(i)  If  x  £  S  is  strictly  e-optimal  for  the  CO(S),  then  x  is  e- efficient  for  the  MOP. 

(ii)  If  x  £  S  is  e-optimal  for  the  CO(5),  then  x  is  weakly  e-efficient  for  the  MOP. 

3.3  Guddat  scalarization 

This  scalarization  developed  by  Guddat  et  al.  (1985)  combines  the  weighted-sum  with  the  constrained- 
objective  scalarization  and  is  defined  as 

m 

G (x°,w):  minimize  ^  Wjfj(x)  subject  to  x  £  S  :=  {x  £  X  :  f(x)  H  f(x°)}, 
i— 1 
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where  x°  E  X  is  a  given  feasible  decision  and  w  E  is  a  given  weighting  parameter. 

Denote  y°  =  f(x°)  and  observe  that  for  this  scalarization  T  =  {y  E  Y  :  y  y°},  U  =  0  and 

n  =  ix  . 

Proposition  3.3  Given  the  MOP  and  e  E  R> ,  let  e  <  wi£i- 

(i)  If  y  E  T  is  strictly  e- optimal  for  the  G(x°,w),  then  y  is  e-Pareto  for  the  MOP. 

(ii)  If  if  ET  is  e-optimal  for  the  G(x°,w),  then  y  is  weak  e-Pareto  for  the  MOP. 

(in)  If  if  E  T  is  e-optimal  for  the  G(x°,w)  with  w  E  then  y  is  e-Pareto  for  the  MOP. 

Proof  In  principle,  each  part  is  proven  as  shown  for  the  weighted-sum  in  Proposition  3.1.  How¬ 

ever,  due  to  the  additional  constraints  in  the  Guddat  scalarization,  we  have  to  verify  that  the 
attainable  outcome  y  E  Y  yielding  the  contradiction  in  each  of  the  three  cases  is  also  attainable  in 
the  Guddat  formulation,  i.e. ,  y  E  T.  We  only  give  the  additional  details  for  (i)  and  then  (ii)  and 

(iii)  follow  analogously. 

Let  if  E  T  be  strictly  e-optimal  for  the  G (x°,w),  and  suppose  that  if  is  not  £-Pareto.  Then  there 
exists  y  E  Y,y  /  y  with  y  >  y  +  e,  yielding 

y<if-e<y<y0 

since  if  E  T  is  attainable  for  G (x°,w).  This  verifies  y  ET,  and  the  proof  is  complete.  □ 

Corollary  3.3  Given  the  MOP  and  e  E  ,  let  e  < 

(i)  If  x  E  S  is  strictly  e-optimal  for  the  G(x°,w),  then  x  is  e- efficient  for  the  MOP. 

(ii)  If  x  E  S  is  e-optimal  for  the  G(x°,w),  then  x  is  weakly  e-efficient,  for  the  MOP. 

(iii)  If  x  E  S  is  e-optimal  for  the  G(x°,w)  with  w  E  then  x  is  e- efficient,  for  the  MOP. 

3.4  Benson  scalarization 

This  scalarization  developed  by  Benson  (1978)  is  defined  as 

m 

B(x°):  minimize  subject  to  f(x°)  +  l  E  f(X)  =  Y.lff  0, 

i=  1 

where  x°  E  X  is  a  given  feasible  decision. 

Observe  that  the  Benson  scalarization  uses  U  =  {l  E  Rm  :  /  =  0}  =  — as  set  of  auxiliary 
variables.  Furthermore,  for  this  scalarization  S  =  {x  E  X  :  f(x)  S  f(x°)},  T  =  {y  eY  :  y  ^  f(x°)} 
and  n  =  X. 
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The  above  formulation  of  the  Benson  scalarization  can  be  rewritten  as 


minimize  yy  fi{x)  —  fi{x°)  subject  to  x  £  S, 

2=1 

where  minxes  {YT=i  /*(x)  “  Mx°)}  =  mina;eS  {YT=  1  fi(x)}  ~  YT=\  fi(x°)-  We  conclude  that 

m 

B(x°)  :  minimize  yy  fi{x)  subject  to  x  £  S  :  =  {x  £  X  :  f(x)  ^  f{x°)} 

i= 1 

is  equivalent  to  the  Benson  problem  B(x°),  and  therefore  that  B(x°)  =  B(x°)  =  G(x°,  1). 

Therefore,  we  can  immediately  derive  that  Proposition  3.3(iii)  and  Corollary  3.3(iii)  also  hold 
true  for  the  Benson  scalarization.  Nevertheless,  we  present  an  alternative  proof  based  on  the  original 
problem  formulation  by  Benson. 

Proposition  3.4  Given  the  MOP  and  e  £  ,  let  e  <  YaL i  £i- 

(i)  If  y  =  y°  +  l  &Y  is  e-optimal  for  the  B(x°),  then  y  is  e-Pareto  for  the  MOP. 

(ii)  If  in  addition  YU i  h  =  0’  then  y°  =  f(x°)  £  Y  is  e-Pareto  for  the  MOP. 

Proof  For  (i),  let  y  =  y°  +  i  £  Y  be  e-optimal  for  the  B(x°)  with  y°  =  f(x°), 

m  m 

li  <  li  +  e  whenever  y°  +  l  £  Y,  l  ^  0, 
i= 1  2= 1 

and  suppose  that  y  is  not  £-nondominated.  Then  there  exists  y  =  y°  +  l  E  Y,  y  ^  y  such  that 
y  >  y  +  £,  yielding 

l  =  y  —  y°<y  —  e  —  y°  =  i  —  effO 

and  thus 

m  m  m  m 

yy  it  >  yy  u + yy  ^  >  yy  h + 

2=1  2=1  2=1  2=1 

in  contradiction  to  the  above.  Hence,  y  is  £-Pareto  for  the  MOP. 

The  second  statement  now  follows  from  the  first  by  observing  that  YJi=i  4  =  0  with  l  ^  0 
implies  that  l  =  0,  and  thus  y  =  y°  =  f(x°).  □ 

Corollary  3.4  Given  the  MOP  and  e  £  M> ,  let  e  <  YJi=\  e*- 

(i)  If  x  £  X  is  e-optimal  for  the  B(x°),  then  x  is  e-efficient,  for  the  MOP. 

(ii)  If  in  addition  YY\  h  =  then  x°  is  e-efficient,  for  the  MOP. 


3.5  Min-max  scalarization 


The  min-max  scalarization,  also  known  as  the  max-ordering  approach  (see  Kouvelis  and  Yu  (1997)), 
is  defined  as 

MM:  minimize  max  {fi(x)}  subject  to  x  £  X. 

Observe  that  for  this  scalarization  S  =  X,  Y  =  T  and  both  U  =  II  =  0. 

Proposition  3.5  Given  the  MOP  and  e  £  R>,  let  e  <  minj=i 

(i)  If  y  €  Y  is  strictly  e-optimal  for  the  MM,  then  y  is  e-Pareto  for  the  MOP. 

(ii)  If  y  €  Y  is  e-optimal  for  the  MM,  then  y  is  weak  e-Pareto  for  the  MOP. 

Proof  We  only  prove  (i)  and  then  (ii)  follows  analogously.  Let  y  £  Y  be  strictly  e-optimal  for 
the  MM, 

max  {yi}  <  max  {y,}  +  e  for  all  y  £  Y,  y  ^  y, 

and  suppose  that  y  is  not  e-Pareto.  Then  there  exists  y  £  Y,  y  /  y  such  that  y  >  y  +  e,  yielding 

max  {yi}  >  max  [ip  +  e,} 

>  max  {yi}  +  min  {ej 

>  max  {y^  +  e 

in  contradiction  to  the  above.  Hence,  y  is  e-Pareto  for  the  MOP.  □ 

Corollary  3.5  Given  the  MOP  and  e  £  Rm,  let  e  <  minj=ii...jm{ei}. 

(i)  If  x  £  X  is  strictly  e-optimal  for  the  MM,  then  x  is  e-efficient,  for  the  MOP. 

(ii)  If  x  £  Y  is  e-optimal  for  the  MM,  then  x  is  weakly  e-efficient  for  the  MOP. 

3.6  Tchebycheff-norm  scalarization 

The  weighted  Tchebycheff-norm  scalarization  of  the  MOP  is  defined  as 

TN(r,  w):  minimize  max  {wiiffx)  —  r*)}  subject  to  x  £  X, 

where  r  £  Mm  is  a  given  reference  or  utopia  point  (see  Steuer  (1986))  and  w  £  R>  is  a  given 
weighting  parameter. 

Observe  that  for  this  scalarization  S  =  X,  T  =  Y,  U  =  0  and  n  =  Mm  x 
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Moreover,  note  that  the  min-max  scalarization  discussed  previously  is  a  special  case  of  the 
weighted  Tchebycheff-norm  formulation.  More  precisely,  if  we  choose  the  reference  point  r  =  0  € 
Mm  and  the  weighting  parameter  w  =  1  G  K>,  then  MM  =  TN(0,1). 

Proposition  3.6  Given  the  MOP  and  e  G  R>,  let  e  <  minj=i)...im{'u;j£j}. 

(%)  If  y  €  Y  is  strictly  e-optimal  for  the  TN(r,w),  then  y  is  e-Pareto  for  the  MOP. 

(ii)  If  y  £  Y  is  e-optimal  for  the  TN(r,w),  then  y  is  weak  e-Pareto  for  the  MOP. 

Proof  We  only  prove  (i)  and  then  (ii)  follows  analogously.  Let  y  G  Y  be  strictly  e-optimal  for 
the  TN(r,  w), 

max  {wiifu  -  ri)}  <  max  {wi(yi  -  r*)}  +  e  for  all  y  G  Y,  y  ±  y, 
and  suppose  that  y  is  not  £-Pareto.  Then  there  exists  y  G  Y,  y  /  y  such  that  y  >  y  +  e,  yielding 

max  {wi(yi  -  ri)}  >  max  {wi(yi  -  r%  +  et)} 

>  max  {Wi{yi  -  n)}  +  min  {w^i} 

>  max  {wi{yi  -  ri)}  +  e 

in  contradiction  to  the  above.  Hence,  y  is  £-Pareto  for  the  MOP.  □ 

Corollary  3.6  Given  the  MOP,  let  e  <  rnirij=i).  ..m{'iLCj£j}. 

(i)  If  x  €  X  is  strictly  e-optimal  for  the  TN(r,w),  then  x  is  e- efficient  for  the  MOP. 

(ii)  If  x  G  X  is  e-optimal  for  the  TN(r,w),  then  x  is  weakly  e- efficient  for  the  MOP. 

3.7  Pascoletti-Serafini  scalarization 

The  direction  method  developed  by  Pascoletti  and  Serafini  (1984)  for  the  MOP  is  defined  as 

PS(r,  v):  minimize  y  subject  to  r  +  yv  —  d  e  f{X)  =  Y,d  E  D, 

where  r  G  Mm  is  a  given  reference  point  and  v  G  D  \  {0}  is  a  given  dominated  direction.  Here 
D  C  Rm  is  a  convex  cone  implying  that  D|DCfl. 

Observe  that  this  scalarization  uses  y  G  M  as  an  auxiliary  variable,  d  G  C  is  a  dominated  slack 
vector  and  n  =  Mm  x  D  \  {0}. 

Proposition  3.7  Given  the  MOP  for  a  convex  cone  D  C  Rm  and  e  G  D,  let  e  G  {5  :  e  —  Sv  G  D} 
and  denote  de  =  e  —  ev  G  D. 
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(%)  If  y  =  r  +  pv  —  d  eF  is  strictly  e-optimal  for  the  PS(r,v),  then  y  is  e-nondominated  for  the 
MOP. 

(ii)  If  y  is  e-optimal  for  the  PS(r,v),  then  y  is  weakly  e-nondominated  for  the  MOP. 

Proof  We  give  the  complete  proof  for  both  (i)  and  (ii).  For  (i),  let  y  with  d  £  D  be  strictly 
e-optimal  for  the  PS(?’,  v), 

p  <  y  +  e  whenever  r  +  yv  —  d  €  Y  for  some  d  G  D, 

and  suppose  that  y  is  not  e-nondominated  for  the  MOP.  Then  there  exists  y  €  Y,y  7^  y  and  d  £  D 
such  that  y  =  y  +  d  +  e,  yielding 


y  =  r  +  pv  —  d  —  d  —  e 
=  r  +  pv  —  d  —  d  —  (ev  +  dt) 

=  r  +  (p  -  e)v  -  (d  +  d  +  de), 

where  d  +  d  +  de  6  D  as  D  is  a  convex  cone.  Now  denote  n  =  p  —  e,  and  then  p  =  p,  +  e  in 
contradiction  to  the  above.  Hence,  y  is  e-nondominated  for  the  MOP. 

For  (ii),  let  y  with  d  £  D  be  e-optimal  for  the  PS(r,  v), 

p<  y  +  e  whenever  r  +  yv  —  d  G  Y  for  some  d  €  D, 

and  suppose  that  y  is  not  weakly  e-nondominated  for  the  MOP.  Then  there  exists  y  €.Y,  y  ^  y,d  € 
int  D  such  that  y  =  y  +  d  +  e.  Since  d  £  int  D,  there  also  exists  7  >  0  such  that  d7  =  d  —  'yv  G  D 
yielding 


y  =  r  +  pv  —  d  —  d  —  e 
=  r  +  pv  —  d  —  ('yv  +  d7)  —  (ev  +  de) 

=  r  +  (p  -  7  -  e)v  -  (d  +  d7  +  de) , 

where  d+d7+de  G  D  as  D  is  a  convex  cone.  Now  denote  y  =  p— 7  — e,  and  then  p  =  y+ 7+e  >  y+e 
in  contradiction  to  the  above.  Hence,  y  is  weakly  e-nondominated  for  the  MOP.  □ 

Corollary  3.7  Given  the  MOP,  D  C  Rm  a  convex  cone  and  e  G  D,  let  e  G  {(5  :  e  —  5v  €  D}. 

(i)  If  x  G  X  is  strictly  e-optimal  for  the  PS(r,v),  then  x  is  e- efficient  for  the  MOP. 

(ii)  If  x  G  X  is  e-optimal  for  the  PS(r,v),  then  x  is  weakly  e- efficient  for  the  MOP. 
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4  Computation  of  e-efficient  solutions 

Based  on  the  results  derived  in  the  previous  section,  we  now  propose  a  generic  procedure  for 
computing  e-efficient  decisions.  The  main  idea  is  that  solving  any  previous  SOP  for  a  (strictly)  e- 
optimal  solution  produces  a  (weakly)  e-efficient  decision  for  the  MOP,  where  the  exact  relationships 
are  summarized  in  Table  1.  Note  that  the  weighted-sum,  Guddat  and  Benson  methods  give  the 
strongest  results  among  the  SOPs  discussed. 

Table  1:  Relationships  between  e-optinrality  and  e-efficiency 


scalarization  method 

strict  e-opt.  =>•  e-eff. 

e-opt.  =>  weak  e-eff. 

e-opt.  =>•  e-eff. 

weighted-sum 

yes 

yes 

yes 

constrained-objective 

yes 

yes 

no 

Guddat 

yes 

yes 

yes 

Benson 

yes 

yes 

yes 

min-max 

yes 

yes 

no 

T  chebycheff-norm 

yes 

yes 

no 

Pascoletti-Serafini 

yes 

yes 

no 

We  derived  that  the  Benson  method  can  be  considered  a  special  case  of  the  Guddat  method. 
Furthermore,  it  is  well  known  that  the  weighted-sum  method  in  general  fails  for  nonconvex  prob¬ 
lems.  Therefore,  we  recommend  that  the  weighted-sum  method  be  applied  for  convex  problems 
and  the  Guddat  scalarization  be  applied  for  nonconvex  problems. 

The  procedure  assumes  that  based  on  prior  knowledge  and  experience,  the  decision  maker  is 
able  to  choose  a  suitable  vector  parameter  e  and  the  scalarization  parameter  n  €  II  required  by  the 
SOP  of  interest.  Clearly,  these  choices  depend  on  the  concrete  problem  to  be  solved,  for  which  more 
specific  advice  can  be  given.  However,  note  that  except  for  the  Pascoletti-Serafini  scalarization, 
each  SOP  implies  an  upper  bound  on  the  set  of  allowable  values  for  e  6  M  such  that  an  e-optimal 
solution  for  the  SOP  is  guaranteed  to  be  (weakly)  e-efficient  for  the  MOP.  For  the  Pascoletti-Serafini 
scalarization,  this  upper  bound  can  be  found  for  the  particular  case  of  D  being  the  Pareto  cone, 
D  =  M>,  and  v  >  0.  Given  that  e  £  {5  :  e  —  Sv  G  M>},  we  obtain 

e  =  max{<5  :  e  —  6v  £  M>  } 

=  max{<5  :  e  ^  <5u} 

=  max{<5  :  et  >  5vi  for  all  *  =  1, ... ,  m } 

=  max{<5  :  £i/vi  >  6  for  allz  =  1, . . . ,  m} 

=  min  {Ei/vi} 
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When  solving  an  SOP  for  e-optimality,  we  recommend  that  the  actual  e  should  be  chosen  as 
the  upper  bound  if  the  largest  tolerance  is  expected.  However,  the  decision  maker  may  choose  any 
other  value  of  e  smaller  than  the  upper  bound  and  hence  gain  control  over  the  e-efficiency  of  feasible 
decisions  for  the  MOP. 

The  proposed  procedure  is  given  in  Table  2. 

Table  2:  Procedure  for  the  computation  of  e-efficient  solutions 

1.  Choice  of  e:  Choose  an  admissible  vector  e  £  Rm. 

2.  Choice  of  the  SOP:  Choose  a  scalarization  method  SOP. 

3.  Choice  of  7r:  Choose  the  scalarization  parameter  it  £  n. 

4.  Computation  of  e:  Compute  the  corresponding  e  £  M. 

5.  Solving  the  SOP:  Solve  the  SOP  for  an  e-optimal  solution. 


Table  3  summarizes  the  required  parameters  n  £  n,  the  admissible  e  £  Rm  and  the  correspond¬ 
ing  upper  bounds  of  e  £  M  for  the  SOPs  examined  in  the  previous  section.  Note  that  in  case  of  the 
min-max  scalarization,  the  third  step  (choice  of  it)  does  not  apply. 


Table  3:  Parameters  for  the  computation  of  e-efficient  solutions 


scalarization  method 

parameter  ir  £  n 

admissible  e 

e  =  e(e,  7r) 

weighted-sum 

w  £  M”1 

£  £  Rm 

e  =  Z2i= i  wi£i 

constrained-objective 

S  £  R™"1 

£  £  R> 

e  =  £i 

Guddat 

(x°,w)  £  X  X 

£  £ 

« =  z2i= i  m£i 

Benson 

x°  £  X 

£  £ 

e  =  £i=i  £i 

min-max 

no  parameter 

£  £ 

e  =  mini=ij...!m{ei} 

T  chebycheff-norm 

(r,  w)  £  Mm  x  R^ 

£  £  Mm 

e  =  mii4=iv..jm{u;j£j} 

Pascoletti-Serafini 

(r,  v)  £  Km  x 

£  £ 

e  =  mini=i  ,...,m{£i/vi} 

5  Example  and  Applications 

We  illustrate  the  procedure  derived  in  the  previous  section  on  the  engineering  example  of  designing 
a  four-bar  plane  truss  as  given  in  Figure  1,  thereby  adopting  the  model  given  in  Stadler  and  Dauer 
(1992)  and  further  discussed  in  Coello  et  al.  (2002). 

The  problem  is  formulated  as  a  biobjective  program  with  the  two  conflicting  objectives  of 
minimizing  both  the  volume  V  of  the  truss  (/i)  and  the  displacement  A  of  the  joint  (/A  subject 
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to  given  physical  restrictions  regarding  the  feasible  cross  sectional  areas  x±,  22,  23,  24  of  the  four 
bars.  The  stress  on  the  truss  is  caused  by  three  forces  of  magnitudes  F  and  2 F  as  depicted  in 
Figure  1.  The  length  L  of  each  bar  and  the  elasticity  constants  E  and  a  of  the  materials  involved 
are  modeled  as  constants. 


Figure  1:  A  four-bar  plane  truss. 

The  mathematical  formulation  of  this  problem  is  given  by 

minimize  j fi(x)  =  L(2x\  +  \[2xi  +  \[2x?,  +  24), 


subject  to  ( F/a )  <  x\  <  3 (F/a), 

V2(F/a)  <  X2  <  3(F/a), 

V2{F/a)  <  23  <  3 (F/a), 

(F/a)  <  24  <  3(F/a). 

Here  the  constant  parameters  are  chosen  as  F  =  10  kN,  E  =  2  x  105  kN/cm2,  L  =  200  cm  and 
a  =  10  kN/cm2.  Since  this  problem  is  convex,  we  chose  the  weighted-sum  method  and  initially 
solved  the  problem  for  five  hundred  weighting  parameters  w  =  (w\  ,1x2),  yielding  five  hundred 
Pareto  outcomes  on  the  Pareto  curve. 

Based  on  the  magnitudes  of  the  found  values  for  truss  volume  (between  1200  cm3  and  2100  cm3) 
and  joint  displacement  (between  0.01  cm  and  0.04  cm),  we  then  allowed  a  tolerance  of  e  1  =  50  cm3 
(200  cm3)  for  the  volume  and  of  £2  =  0.0005  cm  (0.002  cm)  for  the  displacement.  Again  using 
the  weighted-sum  method  with  same  weighting  parameters  w  as  before,  we  solved  the  scalarized 
problem  while  allowing  a  maximal  deviation  of  e  =  w\£\  +  W2£2  from  the  Pareto  optimum  as  a 
stopping  criterion.  Both  the  Pareto  curve  and  the  epsilon-Pareto  outcomes  obtained  in  each  of  the 


2y/2  2y/2  1 


X2 


H - 

23  24 
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two  cases  are  given  in  Figure  2. 


Volume  of  Truss  vs.  Displacement  of  Joint  Volume  of  Truss  vs.  Displacement  of  Joint 


Volume  of  Truss  (cnf)  Volume  of  Truss  (cm) 


Figure  2:  Pareto  and  epsilon-Pareto  outcomes 

As  we  expect,  increasing  the  specified  tolerances  on  truss  volume  and  joint  displacement  from 
comparably  small  to  comparably  large  values  causes  the  e-Pareto  outcomes  to  move  further  from  the 
Pareto  curve.  We  also  observe  that  the  quality  of  e-Pareto  outcomes  varies  significantly  along  the 
Pareto  curve.  The  best  results  are  obtained  for  a  truss  volume  between  1500  cm3  and  2100  cm3  and 
a  joint  displacement  between  0.015  cm  and  0.025  cm.  Moreover,  the  number  of  distinct  outcomes 
generated  is  much  greater  for  the  Pareto  than  for  the  e-Pareto  outcomes. 

While  the  knowledge  of  efficient  decisions  and  Pareto  outcomes  has  proven  to  be  imperative 
for  modern  decision  making,  approximate  efficient  solutions  (or  Pareto  outcomes)  have  not  been 
applied  in  real-life  decision  making.  We  however  believe  that  the  capability  of  computing  approx¬ 
imate  efficient  solutions  may  enhance  and  facilitate  a  decision  making  process  performed  within  a 
multiobjective  programming-based  framework.  We  envision  two  types  of  circumstances  in  which 
approximate  solutions  should  be  relevant. 

First,  robustness  and  sensitivity  of  efficient  solutions  and  Pareto  outcomes  can  be  differently 
approached.  In  the  classical  approach,  an  efficient  decision  x  is  perturbed  by  5x  to  observe  the 
corresponding  change  5y  in  the  outcome  y  and  the  ratio  5x/5y  is  analyzed.  In  the  new  approach, 
perturbation  could  take  place  in  the  objective  rather  than  decision  space.  A  Pareto  outcome  y  could 
be  perturbed  by  Sy  =  e  to  produce  an  e-Pareto  outcome  and  the  related  efficient  and  e-efficient 
decisions  could  be  generated  allowing  thereby  for  the  analysis  of  the  same  ratio  but  in  a  different 
setting. 

The  second  direction  is  related  to  the  relaxation  of  Pareto-optimality  to  e-Pareto  optimality 
so  that  the  set  of  Pareto  outcomes  is  enlarged  with  suboptimal  outcomes.  Such  relaxation  can 
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be  useful  when  dealing  with  a  collection  of  MOPs  for  all  of  which  an  efficient  decision  may  not 
exist.  Here  the  MOPs  yield  a  family  of  Pareto  sets  and  outcomes  that  are  Pareto  for  one  MOP  but 
may  not  be  Pareto  for  another.  This  lack  of  Pareto-optimality  can  be  conveniently  modeled  and 
controlled  with  the  concept  of  e-Pareto  optimality,  and  the  decision  maker  may  gain  control  over 
what  it  takes  for  a  feasible  decision  to  become  e-efficient  for  all  MOPs.  In  this  way,  the  decision 
maker  can  be  more  closely  engaged  in  the  decision  process  not  only  at  the  level  of  choosing  from 
among  efficient  solutions  but  also  at  the  level  of  deciding  solutions’  efficiency. 

6  Conclusion 

Seven  scalarization  methods  are  examined  with  respect  to  their  capability  of  generating  e-efficient 
solutions  for  multiobjective  programming  problems  (MOPs).  It  is  found  that  for  each  method, 
by  appropriate  choices  of  the  involved  parameters,  e-optimal  solutions  for  the  associated  single 
objective  program  become  e-efficient  for  the  MOP.  Quantitative  relationships  between  e  and  e  are 
derived.  An  example  from  engineering  design  illustrates  the  proposed  procedure  for  computing 
e-efffcient  solutions.  Real-life  decision  making  situations  benefiting  from  e-efficiency  are  described 
and  their  resolution  appears  to  be  a  promising  direction  of  further  research. 
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